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Abstract 



We explore how nuclear modifications to the nucleon structure functions, 
shadowing, affect massive gauge boson production in heavy ion collisions at 
different impact parameters. We calculate the dependence of , and 
W~ production on rapidity and impact parameter to next-to-leading order in 
Pb+Pb collisions at 5.5 TeV/nucleon to study quark shadowing at high Q^. 
We also compare our Pb+Pb results to the pp rapidity distributions at 14 
TeV. 
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The heavy ion colhsions at the LHC will be a rich testing ground for hard processes 
which should dominate particle production One of the most promising signatures of 

quark-gluon plasma production at the CERN SPS is J/ip suppression [|,|| which has been 
compared to the Drell-Yan continuum in the lepton pair mass range 2.9 < m < 4.5 GeV 
0. Both J/ip and Drell-Yan pair production are calculable in perturbative QCD. At the 
LHC, quarkonium suppression will be difficult to compare to the dilepton continuum due to 
contributions from cc and bb decays which have large uncertainties in nuclear collisions [Q. 
Since the low mass dilepton continuum is expected to be dominated by bb decays, the 
was suggested as an alternative reference process for quarkonium suppression at the LHC 
II^J^. There are two difficulties with using the Z° as a baseline for quarkonium suppression: 
the large mass differences, mzo ^ mr,mj/^, and the difference in production mechanisms, 
predominantly qq for the and gg for quarkonium. Both these differences are important as 
far as nuclear effects are concerned. However, the differences that reduce the value of the Z° 
as a baseline process are the same that make it an interesting object of study itself — the Z^ 
provides a unique opportunity to study the modifications of the quark distributions in the 
nucleus at high Q^. Therefore, in this paper we examine the possible effects of this shadowing 
on Z° production as well as and W~ production which are also quark dominated. The 
impact parameter dependence of the shadowing effect will also be discussed. 

We further address the issue of how to measure the shadowing effect. Since isospin 
will play an important role in quark-dominated processes, the comparison between Pb-|-Pb 
interactions with and without shadowing is less useful than in gluon-dominated processes 
such as heavy quark production [^. In addition, the first, best, pp data will be at the 
maximum LHC energy of 14 TeV. Therefore we will also present the predicted rapidity 
distributions in Pb+Pb collisions at 5.5 TeV/nucleon, the pp distributions at 14 TeV, and 
the Ph+Ph / pp ratios as a function of rapidity at the two energies. 

The electroweak production and decay channels of the massive vector bosons make them 
excellent candidates for shadowing studies since no hadronic final-state rescattering is pos- 
sible. The Z^ itself, with a 3.37% branching ratio to lepton pairs, will be easily observable 
by reconstructing the peak in the dilepton spectra. Full reconstruction of the leptonic 
decays, l^u, is not possible due to the missing energy given to the undetected neu- 

trino but charged leptons with momenta greater than 40 GeV should be prominent. This 
decay channel has been used at the Tevatron to measure the asymmetry between and 
W~ production since the asymmetry is sensitive to the down to up quark ratio in the proton 
at intermediate values of x and high []TD[ . If the charged leptons from decays can be 



identified in heavy ion collisions, such asymmetry measurements may also be employed at 
the LHC to reduce systematic uncertainties and obtain a more meaningful determination of 
the dependence of quark shadowing in the nucleus. 

The next-to-leading order, NLO, cross section per nucleon for nuclei A and B colliding 
at impact parameter b and producing a vector boson V with mass m at scale Q is 
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where if^^ is proportional to the leading order, LO, partonic ij — > V cross section and 
Q = u,d,s and c. The matrices C" and C"^ contain information on the coupling of the 
various quark flavors to boson V. 

We assume that the parton densities F/^{x,Q^,f, z) can be factorized into x and 
independent nuclear density distributions, position and nuclear-number independent nucleon 
parton densities, and a shadowing function S'^{A, x, Q^, r, z) that describes the modification 
of the nuclear structure functions in position and momentum space. Thus we have 

F^^ix, Q\ r, z) = Pa{s)S\A, x, Q\ f, z)f!'{x, Q') (2) 
Ff{x, Q\ b - f, z') = Pb{s)S^{B, x, Q\ b - f, (x, Q') 

where //^(x, Q^) is the density of parton i in the nucleon and the radial variables s and 



s' are s = \/r^~+~z^ and s' = y\b — r\'^ + z''^. In the absence of nuclear modifications, 
S\A,x,Q'^,r, z) = 1. The nuclear density distribution is described by a Woods-Saxon 
parameterization, 

l + uj{s/RAf 

Pa{s) = po— r, ... . (3) 

1 + exp[(s - RA)/d\ 



Electron scattering data ||rT| are used to fix the parameters Ra, d, u and po- 

Experiments [|12| have shown that the proton and neutron structure functions are mod- 
ified in the nucleus. For momentum fractions x < 0.1 and 0.3 < x < 0.7, a depletion is 
observed in a heavy nucleus relative to a light nucleus such as the deuteron. The low x, shad- 
owing, region and the larger x, EMC, region is bridged by an enhancement at 0.1 < x < 0.3 
called antishadowing. Here we refer to the modification over the entire x range as 'shadow- 
ing' unless otherwise noted. Many theoretical explanations have been proposed, typically 
for only part of the x range such as the very low x or EMC regions. However, as none of 
the models can describe the effect over all x and Q^, we rely on parameterizations of the 
nuclear modifications based on fits to data, as described later. 

Most typical structure function measurements are insensitive to any spatial dependence 
and thus average over the entire nuclear volume. One experiment using a bubble chamber 
found that the structure function does vary spatially but could not determine the dependence 
on impact parameter In a nuclear collision, the impact parameter can be determined 



from the transverse energy production. The influence of the spatial dependence of shad- 
owing on transverse energy production has already been considered |@Jl3- The effects of 
spatially inhomogeneous shadowing on heavy quark 0,|15| , quarkonium, and Drell-Yan PJT6 
production in heavy ion collisions has also been discussed previously. 

We now describe the NLO cross section in Eq. (Q) in more detail. The functions are 
rather simple [jlTj : 

< = f^i(j!,)'^+te«^ (4) 
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where Gp = 1.16639 x 10"^ GeV^, mz = 91.187 GeV, and mw = 80.41 GeV. For Z° 
production by a given flavor i with charge Cj, the sum of the squared vector and axial vector 
couplings is {gyy + {gA)^ = {'^/8){l-4:\ei\xw + 8efxly) where xw = sin^ = l-m^/m|. 



The functions Ajj(x) are universal for all vector bosons, including virtual photons pro- 
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duced in the Drell-Yan process 

cont 

Vg. We have, up to NLO [l^. 



We work in the MS scheme. The NLO correction to 
the qq channel includes the contributions from soft and virtual gluons as well as hard gluons 
from the process qq 



(x) = 6(1 -x) + ^hlS^l _ 4(1 + x)\n( ^) - 8(1 + x) ln(l - x) - 4^^ Inx (6) 
Stt I \myj 1 — X 



Q 



l + x^ 



5(1 



X] 



37r 
6 In 



5)+8C(2) 



m 



16 



1 



In 



Vmf, 



16 



ln(l 



1 — X 



The first delta function is the LO contribution while the NLO contribution is proportional 
to UsiQ"^)- At NLO as(Q^) is calculated to two loops with Uf = 5 active flavors. The last 
three terms are the soft and virtual gluon contributions. The general integral of the 'plus' 
functions in the last two terms is Hi 
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The quark- gluon contribution only appears at 0{as 
At this order Wl 



through the real correction qg qV. 
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For gauge boson production, we take = rriy and all terms proportional to ln(Q^/my) 
drop out. Using the delta functions in Eq. (|l|) we find Xi^2 = {jriy/^/xs) exp(±y). As at LO, 
when Agg(a;) is proportional to 5(1 — x) in Eq. (|^), x'^ 2 = ("^y/\/s) exp(±|/). The rather 
lengthy convolutions of the shadowing functions and parton distribution functions including 
isospin via the proton and neutron numbers are given in the appendix. 

We now define the coupling matrices in Eq. (|lD. The superscripts represent the initial 
(i) and final (f) state quarks or antiquarks while the arguments indicate the orientation of 
the quark line to which the boson is coupled |]T^. The coupling matrices are fairly simple 
for Z° production: C"(gj,gj) = 5ij and C^^{qi,qk) = Sik- With and W~ production, the 
couplings are elements of the CKM matrix. They are nonzero for C^^{qk-iqi) if 6^ + 6^ = ±1 
and for C^^{qk,qi) if = ±1 + e^. In both cases, they take the values |V^j.g;p. Following 
Hamberg et al. [0, we take Vud = cos 9c ~ and Vus = sin 6*0 ~ — with sin 6*0 ~ 0.22. 

We use the MRST HO (central gluon) [0] nucleon parton distributions in the MS scheme, 
shown evaluated at = m\ in Fig. |1|. The valence distributions are somewhat larger than 
the corresponding sea quark distributions at x > 0.1 and extend to higher x values. The 
sea quarks dominate the valence quarks at x ~ 10~^ by a factor of ~ 100. Note also that /| 
is larger than when x > 0.01. The gluon distribution is shown at 1/10 of its magnitude. 
At low X, corresponding to large rapidity, the gluon density is high. 
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The shadowing effect is studied with three parameterizations of the average, homoge- 
neous, shadowing, SKAjXjQ"^) {k = 1 — 3}, measured in nuclear deep-inelastic scattering. 
All the shadowing parameterizations are obtained and evolved at leading order. Since the 
parameterizations are fit to ratios of heavy to light nuclei, the dependence of the param- 
eterizations on both the initial parton densities and the order of the calculation should 
be weak. The first, Si{A,x), assumes that the quark, gluon and antiquark modifications 
are equivalent and includes no evolution |2^. The second, Sl{A,x,Q^), has separate 
modifications for the valence quarks, sea quarks and gluons and includes evolution from 
4 < Q2 < 100 GeV^ [g. The third parameterization, Si{A, x, Q^), is based on the GRV LO 
parton densities. The ratios are evolved over 2.25 < < 10^ GeV^ p3| , p^ assuming 



that Sg"^ = ^3^^ and ^3" = S| while the more massive sea quarks are evolved separately. 
Both the 5*2 and ^3 ratios are evolved to higher using the DGLAP equations pll , p3| , p4| . 
It was shown in Ref. |21[] that including recombination terms in the evolution, as in Ref. ^5 



did not have a large effect on the shadowing ratios, particularly at the x values probed here. 
The initial gluon ratio in 5*3 shows significant antishadowing for 0.1 < a; < 0.3 while the sea 
quark ratios are shadowed. In contrast, 5*2 has less gluon antishadowing and essentially no 
sea quark effect in the same x region. Unfortunately, the evolution of S2 stops below the 
vector boson mass, rendering it less valuable. We show results with all three parameteriza- 
tions because no nuclear DIS data is available at high Q^. Since the S3 parameterization 
includes the most recent nuclear DIS data and is evolved to scales compatible with the vector 
boson masses, the 5*3 results should perhaps be favored. 

The shadowing ratios in a lead nucleus compared to a proton are shown in Fig. 0. The 
effects of shadowing on the valence quarks is strongest with 5*1 since all quarks are treated 
equally. The S2 and S3 valence ratios are rather similar in magnitude although the anti- 
shadowing range is broadest for 5*3, 0.01 < a; < 0.3 and the 5*3 ratio is lower than the 5*2 ratio 
at low X. The 5*1 and 5*2 sea quark ratios are very similar when x < 0.1. Then the 5*2 ratio is 
essentially unity until x > 0.3. The 5*3 ratios are all larger than the 5*1 and 5*2 ratios, even at 
small X, due to evolution. It is most interesting to note the differen_ce between the light and 
strange sea ratios in the S3 parameterization. The ratios S3 and S3 show no antishadowing 
effect but instead decrease when x > 0.1 while S^ and S^ are typically larger over all x and 
are antishadowed when 0.01 < a; < 0.2. The antishadowing in the charm distribution is 
larger even than for the gluons. Since all 5*3 sea ratios are equivalent at = 2.25 GeV^, 
the difference is solely the effect of evolution. Note also that the gluon shadowing ratios 
are typically larger than the large sea quark ratios over all x. At = m| the strong 
antishadowing in S3 has essentially disappeared and is no larger than that of S2 although 
the antishadowing region of S3 is broader, from 0.005 < x < 0.2. 

Nuclear shadowing should depend on spatial position of the partons in the nucleus as well 
as on their momentum. Most models predict some form of spatial dependence, according 
to the origin of the shadowing effect. Typically, the spatial dependence can be expected to 
take two forms, either proportional to the local nuclear density, Eq. (^, or the path length 
of the parton through the nucleus. Both will be discussed below. 

When the parton density is high, partons in one nucleon can interact with those in 
neighboring nucleons, recombining to lower the parton density |^|. In this case shadowing 
is proportional to the local nuclear density 0,113]. Then 
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^^ws = Sl{Ax,Q\r,z) = l + Ny,s[Sl{A,x,Q') - 1]^ , (9) 

Po 

where iVws is a normalization constant chosen so that (l/A) J d'^spA{s)SlT^s = SI. At large 
radii, s ^ Ra, the medium modifications weaken and the nucleons behave as though they 
were free. At the center of the nucleus, the modifications are larger than the average value 
determined from nuclear DIS. 

It has also been suggested that shadowing stems from multiple interactions of the incident 



parton [^. In this picture, parton-parton interactions are longitudinally distributed over 



the coherence length, = l/2Tnj\ix, where m^v is the nucleon mass [Q. When x < 0.016, 
Ic > Ra for all nuclei and the interaction of the initial parton is delocalized over the entire 
nuclear path, thus interacting coherently with all target partons along the distance Ic- For 
small X, shadowing depends on the longitudinally-integrated nuclear density at transverse 
distance r and the spatial dependence can then be parameterized as 



Sl^iA, X, Q\ r, z) = l + N,{Sl{A, x, Q') - (10) 

where the normalization is again defined by (l/A) J d?rdzpA{s)S\p = SI with Np > A^ws- 
However, at large x, Ic -C Ra and shadowing is again proportional to the local density so 
that Eq. (^) corresponds to the large x limit of the multiple scattering formulation. 

There are some problems with implementing the multiple scattering picture in nuclear 
collisions. While traversing the formation length, both the initial- and final-state par- 
tons may undergo multiple interactions, reducing the effective Ic, similar to the Landau- 
Pomeranchuk-Migdal effect [^]. In addition, the idealized picture of a single initial parton 
incident on a static nucleus is inappropriate in heavy ion collisions since many interactions 
occur simultaneously, increasing the density in the path of the initial parton. A cascade 
approach cannot resolve the difficulty because non-local depictions of these collisions are 
Lorentz frame dependent |3^. Finally, since the parton densities are distributed over an 



x-dependent distance, baryon number is not locally conserved even if the valence quarks are 
considered to be fixed spatially. Given the difficulties with the multiple interaction picture 
as well as those of matching the two spatial dependencies according to Ic at each x, we only 
present specific results for the local density model, Eq. (|). 



Other mechanisms of shadowing effects in the EMC region such as nuclear binding 31 



and rescaling |]32|j33[| have also been suggested but can explain only part of the observed 
effect |^4|. We note that these models would also predict some spatial dependence. 

We first show that our shadowing results do not depend strongly on the order of the 
calculation. The ratio of the NLO to LO cross sections, both calculated with the MRST 
HO distributions, is often referred to as the K factor. The K factor is given as a function 
of rapidity with no shadowing for Pb+Pb collisions at 5.5 TeV in Fig. ^ for all three vector 
bosons. The K factor is ~ 1.13 up to ?/ = 3. It grows larger as the edge of phase space 
is approached since large, positive, y corresponds to low X2 where the gluon density is high 
and the qg channel becomes more important. The qg channel contributes ~ 15% of the 
total vector boson cross section, with or without shadowing. The K factor increases faster 
with y for the Z° because the higher mass means that the phase space for production 
is exhausted at lower rapidities than production, leading to larger qg contributions at 
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high yz- The K factors are quite similar when shadowing is included and differ from those 
without shadowing by ~ 1% at ?/ = 0. At high rapidities, a larger effect might be expected 
from the gluons but, as seen in Fig. 0(c), the effects of shadowing on the gluon distributions 
are not as strong as those of the sea quarks at large Q^. Thus the difference in the K 
factors between the calculations with and without shadowing is only marginally larger at 
high rapidity, up to ~ 4% at ?/ = 4. 

In Fig. ^, we show the ratio of production in Pb+Pb collisions with and without 
shadowing at both LO and NLO. The results are independent of the order of the calculation, 
even at large rapidities. This is not surprising since we have shown that the differences 
between shadowing ratios at LO and NLO are trivial for virtual photon production via the 
Drell-Yan process 0, even for pairs with m < mx. At the higher scale of and 
production, the approximation should be even better because the K factor is smaller. 

We now calculate the NLO and W~ cross sections in nuclear collisions. Table | 

gives the total cross sections in the CMS and ALICE central acceptances, \y\ < 2.4 and \y\ < 
1 respectively, at the LHC. The cross sections are larger than the virtual photon mediated 
Drell-Yan cross sections at lower masses 0. The results, given for Pb+Pb collisions, are 
integrated over impact parameter and presented in units of nb/nucleon pair. We note that 
with the normalization of S'ws? the impact-parameter integrated cross section is unchanged 
when the spatial dependence is included. The next-to- next-to-leading order, NNLO, 
and Z*^ total cross sections have also been calculated fl^ . The K factors obtained from the 
ratio of the 0{aa1) to 0{a) cross sections differs by ~ 1% from the 0{aas) to 0{a) K 
factor shown in Fig. |^. Thus changes in the total cross sections between NLO and NNLO are 
on the few percent level even though the vector bosons can be produced in the gg channel 
as well at NNLO, because asimy) ~ 0.116. The effects on the shadowing ratios should be 
even smaller, see Fig. ^. 

We have checked how our results depend on the chosen set of parton densities. Using 
the CTEQ5M densities the total cross sections in Table | increase by ~ 5% over the 
MRST HO results. However, the K factors and shadowing ratios in Figs. |^ and ^ change 
by less than 1%. Thus our shadowing results are essentially independent of parton density. 

In Table we show the expected rate in nucleus-nucleus collisions at 6 = 0, N{S = S^) = 
(TAr7v7pbPb(0)i^PbPb<^('S' = 5*^) with i^pbpb = 1/nb in a one month (10^ s) LHC run, (Ttvat = 60 
mb at LHC energies, and TpbPb(O) = 30.4/mb. The absolute numbers in the experimental 
acceptances are large but do not reflect the measurable decay channels. Including the 3.37% 
lepton pair branching for Z^ decays reduces the number produced with no shadowing, 5* = 1, 
to 990 in CMS and 425 in ALICE. The 10% lepton branching ratio for W+ and W~ leaves 
nearly 4600 observable decays in CMS and 1980 in ALICE. 

Figures ^ and ^ compare the ratios of Z° and production in Pb+Pb collisions with 
the three shadowing parameterizations to Pb+Pb collisions with no shadowing as a function 
of rapidity. The isospin effects wash out the differences between the and W~ distri- 
butions in the ratios so that the results are essentially identical for the two charged vector 
bosons. Therefore the ratios are shown only for the W~^. The results are given for several 
impact parameter bins, the most central bin, b < 0.2Ra, an intermediate impact parameter 
bin around b ~ Ra, and a peripheral bin around b ~ 2Ra- It is clear that by neglecting the 
impact parameter dependence of shadowing, one may overestimate the effect in peripheral 
collisions, an important point if using the Z^ as a baseline in different transverse energy 
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bins. The integration over all impact parameters is equivalent to the average shadowing, as 
expected from the normalization of Eq. (P). Although the results are shown using the local 
density approximation, the parameterization of the spatial dependence for a long coherence 
length, Eq. (p!oD, differs only marginally. In central collisions, the difference between the two 
parameterizations is less than 1% while in the most peripheral bin, it is 3 — 6%. The largest 
differences occur in regions with the strongest shadowing modifications. Thus the calcu- 
lations are rather insensitive to the exact spatial parameterization, suggesting that heavy 
ion collisions cannot distinguish between different dependencies, only between homogeneous 
and inhomogeneous shadowing. 

The ratios are rather similar for all vector bosons. The 5*1 and 5*2 ratios are approximately 
equal as a function of rapidity, presumably because the evolution of the 5*2 parameter- 
ization ends at = 100 GeV^. The calculations cover the entire rapidity range of vector 
boson production. At yz ~ 0, = X2 = 0.017, in the low x region. As rapidity increases, 
Xi increases, going through the antishadowing region and the EMC region with xi ~ 0.33 
at yz = 3. The kink in the Si ratio at ?/ ~ 2.2 is an artifact resulting from the rather 
sharp transition between the shadowing and EMC regions at a; ~ 0.15, see Fig. 0. Note that 
the larger rapidity coverage of CMS makes the EMC region accessible in this measurement. 
When yz — > 4, a;i — > 1, entering the "Fermi motion" region and causing the upturn of the 
ratios at large yz- Note also that at large Xi, the valence quarks dominate. While increasing 
yz (xi) traces out the large x portion of the shadowing curve, the low x part of the shadow- 
ing regime is accessible from X2 with growing yz- At = 3, X2 ~ 8 x 10^^, in a range where 
shadowing saturates in Si and 5*2. There is no saturation built into the 5*3 parameterization, 
causing a steeper decrease in the ratios for large yz with this parameterization than with Si 
and 5*2. In addition, the 5*3 sea quark shadowing is never as strong at low x as for Si and 
5*2 so that these two parameterizations are both more strongly shadowed overall. The Z° 
ratios are all slightly higher than those for because the larger mass of the Z° results in 
xz ~ l.lxw- 

We also point out that the large vector boson masses do not allow us to restrict ourselves 
only to an x region where the coherence length is always larger than so that the multiple 
interaction approach could be used without having to match the spatial dependence across 
X boundaries. While the target parton is at relatively low x, the projectile parton, also 
affected by shadowing, is at relatively high x where Ic is small. 

The shadowing ratios are fairly simply traced out for vector boson production, especially 
at leading order since the fixed boson mass defines x at any y whereas Drell-Yan shadowing 
effects are smeared over the mass interval. However, the ratios shown in Figs. |^ and ^will not 
be accessible experimentally due to the nuclear isospin. The comparison must be made to pp 
interactions, preferably at the same energy to retain the same x values. This ideal situation 
may not be realized for some time at the LHC. Therefore in Figs. we show the Pb+Pb 
rapidity distributions with and without homogeneous shadowing as well as the distributions 
from pp collisions at 14 TeV for all three vector bosons. The Pb-|-Pb cross sections are given 
per nucleon pair for a more direct comparison. The higher energy extends the available 
vector boson rapidity space by one unit. The distributions in Fig. |^ have a plateau over 
two units of rapidity. The pp distribution in Fig. |^ rises over the first several units of 
rapidity, followed by a decrease as the edge of phase space is approached. This is due to the 
increasing importance of valence quarks at large y (large Xi). The effect appears for 
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production in pp collisions because the u valence quarks carry more momentum than the d 
valence quarks, see Fig. |l] (a). On the other hand, the W~ pp distribution always decreases 
with rapidity. In Pb+Pb collisions, there is a slight increase in W~ production with rapidity 
instead of a decrease while the distributions are either flat or decreasing. This increase 
in W~ production shown in Fig. ^ is due to the neutron excess in Pb+Pb where, in nn 
collisions, W~ production proceeds dominantly through f^f^i^ fjfu)- Likewise, the rise in 
production in Pb+Pb relative to pp collisions disappears because of the neutron content 
of the nucleus. 

Finally, the ratios of the Pb+Pb/nucleon pair to pp cross sections are shown in Figs. 
|12| for homogeneous shadowing. Due to the higher pp cross sections, the ratios are lower than 
those to = 1 at 5.5 TeV shown in Figs. ^ and p. Since both the Pb+Pb and pp rapidity 
distributions are rather fiat, the ratios in Fig. ^ are also fiat to ?/ ~ 1.5. The rise in the 
pp distribution shown in Fig. |^ causes the Pb+Pb to pp ratios to decrease with rapidity 
over all y in Fig. 11. However, the increase in W~ production due to neutrons in Pb+Pb 
collisions, barely visible in Fig. ^ is apparent in the W~ Pb+Pb to pp ratio in Fig. |T2|. It 
should still be possible to distinguish between the shadowing parameterizations and study 
quark shadowing at = my, particularly since the 14 TeV pp data will be available 
with higher statistics. Note that comparing the pp results to the Pb+Pb calculations with 
inhomogeneous shadowing would result in slightly lower ratios in central collisions and higher 
ratios in peripheral collisions, as expected from Figs. H and 

Once the basic nuclear shadowing effects on vector boson production have been under- 
stood, they can perhaps be used to study other medium effects in heavy ion collisions by 
comparing the leptonic and hadronic decay channels. The hadronic decays of the vector 
bosons, ~ 70% of all decays of each boson, may be more difficult to interpret. While the 
width of the Z° decay to is not expected to be modified in the quark-gluon plasma due 
to the weak coupling ||3^, the has a 2.49 GeV total width and will decay in any quark- 
gluon plasma to two jets through Z° — »• — ^ jet + jet in ~ 0.1 fm. Therefore, the decay 
jets could be modified in the medium which may still be progressing toward thermalization 
and will be subject to rescattering and jet quenching. Thus a comparison of a reconstructed 
Z^ in the dilepton channel where no nuclear effects are expected and medium-modified jets 



should result in a broader width in the qq channel than the rl~ channel In addition, 
the Z° could be used to tag jets through the qq — > Z'^g and gq Z^q channels to study jet 
properties in the quark-gluon plasma [If. 
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Appendix 



In AB collisions, the cross section per nucleon must include the nuclear isospin since, in 
general, ^ o"^ ^ ^ cr^„. We give the convolution of the nuclear parton densities in 
Eq. g), including only the couplings. We take = f^^, f^^ = f^, = /|, and = /|. 
All other distributions are assumed to be identical for protons and neutrons. The proton 
and neutron numbers in nucleus A are Za and Na- To be concise, we define 

S'^{A,x)mx,Q') = S--{A,x)fl^ + 5^(A,a;)/| (Al) 

S\A,x)ax,Q') = S''-{A,x)fl + S\A,x)f^^ (A2) 
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where we have abbreviated the shadowing functions as S^{A,x). 
We begin with the qq channel For Z° production, we have 

^_5^(^,x05^(s,x2)/,^(xi,g^)/|^(x2,g^)c^'(?„?,)[W' + i9\r] (A3) 

i,jeQQ 

1 8 '-{2 

= -[1 - -xw + -xl^] {S-{A,xr)S''{B,X2) Q^) + 7V^/:(xi, Q^)} 

X {ZBfiix^, Q') + NBmx2, Q')} + 2ABS''{A, x,)S%B, a;^)/^^!, Q')/I(a;2, Q')) 

+ ^[1 - + ^x^v^] {S''iA,x,)S'{B,X2) {ZAf!{x^,Q') + NAfS{x^,Q')} 

X g^) + NBf^{x2, Q')} + 2ABS^{A, x,)S^{B, a;^)/^ (xi, g^)/|(a;2, Q^)) 

+ [xi ^ X2, A ^ B] . 

Note that for production, we have also included the square of the vector and axial vector 
couplings since these depend on the quark charges. The qq convolution is 

Y^j\A,x,)S^{B,X2)f,';{xr,Q')f-^^{x2,Q')C'%,q^^ (A4) 

= cos^ 9c {S-{A, x,)S\B, X2) {ZaFu{x,, Q') + Na^x,, Q')} 

X {ZBflix2, Q') + NBf^ix2, Q') } + ABS^iA, x,)S%B, Q')f!ix2, Q')) 

+ sinHc {S^iA,x^)S^iB,X2) {ZAmx,,Q') + NAf:ix,,Q')} BfUx2,Q') 

+S'iA, x,)S%B, X2) {ZAf^ixi, Q') + NAf^ixu Q')}Bf^ix2, Q')) } 
+ [xi ^ X2, A ^ B] . 

Finally, the W~ qq convolution is 

J2_S\A, x,)S^{B, X2)f^{xu Q')f£ {X2, Q')C'%, qj) (A5) 

= cos^ Oc {S'^iA, x,)S''{B, X2) {ZAfU^i, Q') + Na^x,, Q') } 

X {ZBfAx2, Q') + NBnix2, Q')} + ABS^iA, x,)S\B, a;2)/f (xi, Q^)fl{x2, Q^)) 

+ sin^ dc x,)S\B, X2) {ZaM{x,, Q') + Na^x^, Q')}Bf^ix2, Q') 

+S''iA, x^)S~%B, X2) {ZAfAxu Q') + A^a/,"(xi, Q')} 5/1(^2, Q')) } 
+ [xi '(-^ X2, A B] . 

We now turn to the qg channel. The convolution for Z^ production is 
E_ (S\A, x,)S%B, X2)f^^{x„ Q')f^{x2, Q') 

i,keQQ 

+ [xi ^X2,A^ B]y'\q„ q^Wvf + (g\f] (A6) 
= BS<^{B,X2)f^{x2,Q')[l[l - ^xw + yx^] (S'^iAx,) {ZAf^xuQ') + NAf:{x,,Q')} 
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+ ^[l-^xw + Ixl,] (S'^iA, X,) {ZAfAx,, Q') + NAfSix,, Q')} 

+ S'(A,x,) {ZBfl{x2,Q') + NBf^{x2,Q')} + 2AS%A,x^)f^{x„Q')) | 
+ [xi X2, A B] . 

For production in the qg channel, we have 



J2_[^S\A,x,)S%B,X2)f,';{x„Q')f^{x2,Q') + [x, ^X2,A^ B]jC^'{q„qk) (A7) 

i,keQQ 

= BS^{B,X2)f^{x2,Q')[s'^{A,x,) {ZAmx,,Q') + NAf:{x,,Q')} 

+S'{A,x^) {ZBf^{x2,Q') + NBf^ix2,Q')} + A{s^{A,x^)f|{x^,Q') + S^iA,x,)f^{x^,Q')} 
+ [xi X2, A B] . 

Now the couphngs do not enter exphcitly for and production because each distri- 
bution is multiphed by ( cos^ 6'c + sin^ 6'c)- Finally, the qg convolution for W~ production 
is 

^_ i^S%A,x,)S^{B,X2)f^';{x„Q')f^{x2,Q') + [x, ^X2,A^ Bj^C'^q^q,) (A8) 
= BS%B, X2)f^{x2, Q') [s'^iA, X,) {ZAfiix^, Q') + TV^/^ (xi, Q')} 

+S\A, X,) {ZBfAx2, Q') + NBfS{x2, Q')} + A {S'{A, x^/f (xi, Q') + S'^iA, x,)fl{xu Q')} 
+ [xi ^ X2, A ^ B] . 
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TABLES 



Detector a{S = 1) (nb) a{S = Si) (nb) a{S = S2) (nb) a{S = S3) (nb) 

CMS 16.10 11.37 11.22 14.92 

ALICE 6.93 4.87 4.93 6.56 

CMS 25.18 17.39 17.08 23.23 

ALICE 10.84 7.39 7.45 10.19 

w- 

CMS 26.63 18.39 18.12 24.58 

ALICE 11.21 7.64 7.73 10.54 

TABLE I. Vector boson production cross sections in units of nb per nucleon pair in Pb+Pb 
collisions at 5.5 TeV/nucleon calculated with the MRST HO parton densities. Full azimuthal 
coverage is assumed. The corresponding pp cross sections at 14 TeV are a^" = 35.44 nb (CMS), 
14.94 nb (ALICE), a^^ = 60.50 nb (CMS), 24.76 nb (ALICE) and = 52.95 nb (CMS), 22.88 
nb (ALICE). 



Detector 


N{S = 


= 1) 


N{S = Si) 


N{S = 


S2) 


N{S = 


S3) 


CMS 


2.94 X 


lO'' 


2.07 X 10^ 


2.05 X 


10^ 


2.72 X 


10" 


ALICE 


1.26 X 


lO'' 


8.88 X 10^ 


8.99 X 


10^ 


1.20 X 


10" 








W+ 










CMS 


4.59 X 


10^ 


3.17 X 10'' 


3.12 X 


10^ 


4.24 X 


10^ 


ALICE 


1.98 X 


lO"* 


1.35 X 10^ 


1.36 X 


10^ 


1.86 X 


10" 








w- 










CMS 


4.86 X 


lO'' 


3.35 X 10^ 


3.31 X 


10^ 


4.48 X 


10" 


ALICE 


2.04 X 


lO'' 


1.39 X 10^ 


1.41 X 


10^ 


1.92 X 


10" 



TABLE II. Number of vector bosons produced at 6 = in a one month (10^ s) Pb+Pb LHC 
run at 5.5 TeV/nucleon. Note that no decay branching ratios have been included. 



13 



FIGURES 




X X 

FIG. 1. The MRST HO proton parton distribution functions evaluated at Q = mz- The up 

(sohd) and down (dashed) valence distributions are given in (a) while the up (lower solid), down 
(dashed), strange (dot-dashed) and charm (dotted) sea quark distributions are shown in (b), along 
with the gluon distribution (upper solid), reduced by a factor of 10 for comparison. 
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FIG. 2. The homogeneous shadowing parameterizations used in our calculations, evaluated at 

q2 _ rv?2. Valence shadowing is shown in (a) for the 5*1 (solid), 5*2^ (dashed), and 6*3 (dot-dashed) 
parameterizations. Sea quark shadowing is shown in (b) for S\ (solid), (dashed), 5*3 = 5*3 
(dot-dashed), (dotted) and (dot-dot-dot-dashed). Gluon shadowing is shown in (c) for 5*1 
(solid), 51 (dashed) and (dot-dashed). 
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FIG. 4. The shadowing results at LO and NLO are compared. The NLO results are given in 
the dashed, Si, dot-dashed, ^2, and dotted, 53, lines. The LO shadowing ratios for Si, circles, S2, 
squares, and 53, diamonds, are also shown. 
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FIG. 5. The rapidity distributions, relative to 5' = 1 for Pb+Pb collisions at the LHC, 
calculated with the MRST HO distributions. Central, h < 0.2i?^, semi-central, 0.9i?A <h < \.\Ra-, 
and peripheral, 1.%Ra < b < 2.1Ra impact parameters are shown along with the integral over all b. 

The homogeneous shadowing results are given in the dashed, Si, dot-dashed, 52, and dotted, ^3, 
lines. The inhomogeneous shadowing ratios for wSj circles, 52WS) squares, and S'swSj diamonds, 
are also shown. 
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FIG. 6. The rapidity distributions, relative to 5 = 1 for Pb+Pb collisions at the LHC, 
calculated with the MRST HO distributions. Central, h < 0.2Ra, semi-central, O.QRa < b < I.IRa, 
and peripheral, 1.9Ra <b < 2.1Ra impact parameters are shown along with the integral over all b. 

The homogeneous shadowing results are given in the dashed, ^i, dot-dashed, 5*2; and dotted, ^3, 
lines. The inhomogeneous shadowing ratios for 5iwS) circles, 52 wS) squares, and 5*3 wSj diamonds, 
are also shown. 
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FIG. 7. The Z° rapidity distributions in pp and Pb+Pb collisions, calculated with the MRST 
HO distributions. The upper solid curve is the pp result at 14 TeV while the lower solid curve is 
the Pb+Pb distribution at 5.5 TeV/nucleon pair with no shadowing. The homogeneous shadowing 
results for Pb+Pb collisions are given in the dashed, Si, dot-dashed, S2, and dotted, Ss, lines. 
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FIG. 8. The rapidity distributions in pp and Pb+Pb collisions, calculated with the MRST 
HO distributions. The upper solid curve is the pp result at 14 TeV while the lower solid curve is 
the Pb+Pb distribution at 5.5 TeV/nucleon pair with no shadowing. The homogeneous shadowing 
results for Pb+Pb collisions are given in the dashed, Si, dot-dashed, S2, and dotted, Ss, lines. 
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FIG. 9. The W rapidity distributions in pp and Pb+Pb collisions, calculated with the MRST 
HO distributions. The upper solid curve is the pp result at 14 TeV while the lower solid curve is 
the Pb+Pb distribution at 5.5 TeV/nucleon pair with no shadowing. The homogeneous shadowing 
results for Pb+Pb collisions are given in the dashed, Si, dot-dashed, S2, and dotted, Ss, lines. 
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FIG. 10. The ratios of the rapidity distributions in Pb+Pb cohisions relative to pp colhsions, 
calculated with the MRST HO distributions. The solid curve is the ratio without shadowing. The 
homogeneous shadowing results are given in the dashed, ^i, dot-dashed, S2, and dotted, ^3, lines. 
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FIG. 11. The ratios of the rapidity distributions in Pb+Pb coUisions relative to pp coUi- 
sions, calculated with the MRST HO distributions. The solid curve is the ratio without shadowing. 
The homogeneous shadowing results are given in the dashed, ^i, dot-dashed, S2, and dotted, S3, 
lines. 
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FIG. 12. The ratios of the W~ rapidity distributions in Pb+Pb coUisions relative to pp coUi- 
sions, calculated with the MRST HO distributions. The solid curve is the ratio without shadowing. 
The homogeneous shadowing results are given in the dashed, ^i, dot-dashed, S2, and dotted, S3, 
lines. 
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